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Abstract Although only a small fraction of stars end their lives as supernovae, 
all supernovae leave behind a supernova remnant (SNR), an expanding shock wave 
that interacts with the surrounding medium, heating the gas and seeding the cosmos 
with elements forged in the progenitor In this chapter, we introduce the basic prop- 
erties of galactic and extragalactic SNRs ($2). We summarize how SNRs evolve 
throughout their life cycles over the course of ~ 10° yrs (83). We discuss the var- 
ious morphological types of SNRs and discuss the emission processes at various 
wavelengths.($4). 
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2 Introduction 


Supernova (SN) explosions are among the most energetic events in the universe 
since the Big Bang, releasing more energy (~ 10?! — 10°? ergs) than the Sun will 
release over its entire lifetime. They are the cataclysmic ends of certain types of 
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stars, and are responsible for seeding the universe with the material necessary to 
form other stars, planets, and life itself. We owe our very existence to generations of 
stars that lived and died billions of years ago, before the formation of our Sun and 
solar system. The processes of stellar evolution continue to occur today, with typical 
galaxies like the Milky Way hosting a few supernovae per century, on average. The 
study of supernovae and the role they play in shaping the evolution of star systems 
and galaxies is truly an exploration of our own origins. Supernova remnants (SNRs), 
the expanding clouds of material that remain after the explosion, spread elements 
over volumes of thousands of cubic light-years, and heat the interstellar medium 
through fast shock waves generated by the ejecta from the star. 

There are currently ~ 300 SNRs cataloged in our Galaxy [1], although radio 
surface brightness studies predict that only half of SNRs have been identified [2]. 
Recent X-ray observations with good spatial resolution and large effective area al- 
low us to detect many SNR samples in nearby galaxies as well, such as Magellanic 
Clouds[3], M31[4], and M33[5]. 


3 Evolution of supernova remnants 


The general picture of an SNR is that of a cloud of expanding material ejected from 
the star, rich in heavy elements like O, Si, and Fe, that drive a shock wave into the 
interstellar medium (ISM). This shock waves heat the interstellar gas it encounters 
to millions of degrees, creating a highly-ionized plasma. Figure 1 shows Cassiopeia 
A in X-rays, an example of a young (~ 350 yr) SNR in our own galaxy. Although 
Cassiopeia A is known to have resulted from a core-collapse (CC) SN (see Chapter 
on SN), it is generally difficult to tell the type of SN only by looking at the remnant. 
SNRs remain visible for thousands, often tens or hundreds of thousands of years 
before dissipating their energy into the ISM. The life of a SNR can be thought of as 
consisting of four phases. 


3.1 Free expansion phase 


The ejected gas mainly determines the evolution of SNRs in the initial phase. Since 
the pressure of surrounding matter is negligible in comparison, the exploded ejecta 
expands effectively without deceleration. As a result, this phase is called the “free 
expansion" phase. 

The characteristic (or typical) speed of the ejecta expansion, v,, is given by 


1) 
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Fig. 1 Cassiopeia A, a young remnant in our Galaxy, as seen at X-ray energies with Chandra. 
Credit: NASA/CXC/SAO 


where E and M,; are the kinetic energy of the explosion and the ejecta mass, respec- 
tively. With fiducial values of Egin = 1 x 10?! erg and Me jį = 1Mo, a typical velocity 
of the ejecta is ~ 10^ km s~!. While this velocity for the ejecta is quite high, it 
is also a bulk velocity with very little randomization in the velocity of the various 
fluid elements without deceleration (This is the difinition of “free expansion"). 
Thus, there is no heating and the ejecta temperature is actually quite cold in the 
beginning, and the remnant is generally faint at all wavelengths. The exception is 
the radiation from nonthermal synchrotron emission from high energy electrons 
accelerated in the shock (see Chapter *acceleration"). 

Virtually all SNRs emit in radio waves, where the emission is entirely nonthermal 
in origin, resulting from synchrotron radiation from relativistic electrons spiraling 
around magnetic fields in the shock. Synchrotron emission is characterized by a 
featureless power-law spectrum, where the radio flux, Sy is given by Sy œ v ^, 
where v is the frequency and a is the spectral index, which depends on the energy 
distribution of the electron population. For more detail, see Chapter “‘accelera- 
tion". In Figure 2, we show an example of a remnant, G1.94-0.3, which is still in 
the free expansion phase [6]. The acceleration of electrons up to GeV energies is 
sufficient to produce radio synchrotron emission. 

A small population of SNRs are apparently able to accelerate electrons to much 
higher energies (> 10 TeV), sufficient to produce synchrotron radiation at X-ray 
energies [7, 8]. G1.9+0.3 also falls into this category, and Figure 2 also shows an 
image of the remnant that is essentially entirely from nonthermal synchrotron emis- 
sion. For more details regarding the particle acceleration mechanisms, see the review 
by [9]. 
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Fig. 2 Left: SNR G1.9+0.3, the youngest remnant in the Galaxy, seen in radio waves from the 
VLA in 1985. Right: The remnant, in 2007, seen in X-rays with Chandra. In the 22 years between 
these images, the remnant expanded substantially, allowing the time since the original supernova 
explosion (about 140 years) to be estimated. Credit: X-ray (NASA/CXC/NCSU/S.Reynolds et al.); 
Radio (NSF/NRAO/VLA/Cambridge/D.Green et al.). 


3.2 Adiabatic expansion phase 


At the end of the free expansion phase, the expansion begins to decelerate due to 
the forward shock sweeping up a non-negligible amount of material in the ISM. An 
estimate for when this transition occurs can be obtained by calculating the radius, 
R$, at which the swept-up mass equals the ejecta mass: 


4 3 
— Re po = Me; , (2) 


3 
Ma; NP n -1/3 
T ij aw 

i s(a) Gr cm? ) (po) , 6) 


where po and n are mass and number densities in the cgs unit. Assuming typical 
values of supernova explosion and ISM, it generally takes a few hundred years for 
the radius of the SNR to reach R? (though this can vary by orders of magnitude if 
the SN explodes into particularly low or high density regions of the ISM). 

This phase of deceleration of the forward shock wave also leads to the formation 
of a "reverse shock," which propagates back into the ejecta towards the center of the 
SNR. Similarity solutions for the reverse shock are discussed in Hamilton & Sarazin 
[10, 11]. 

In this phase, the evolution of the SNR radius as a function of time can be written 
as 

R(t)ect" , (4) 
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where m is the expansion parameter defined by [12]. At the end of the free expansion 
phase, m ~ 1, by definition, but during this second phase, m gradually decreases to 
m — 0.4 [13]. 

When the SNR shock starts decelerating, radiative cooling is still negligible and 
the expansion is adiabatic. For the most simple case of a uniform ISM, the evolution 
can be written as the results of the similarity solution by [13] (thus, this phase is 
often referred to as the “Sedov phase"). 

Assuming that a large amount of energy is instantaneously released in a small 
volume, expansion can be characterized by only two parameters, po (the mass den- 
sity of the ISM) and Ekin (the kinetic energy of the explosion). We can define a 
non-dimensional parameter Ó as 


Er? —1/5 
=R; FA , 
S ( po ) à 


where R, and t represent the radius and age of the SNR. R, can be written as 


E 2 1/5 
[e (=) , (6) 
po 


where Co is determined from the energy conservation equation. Using Co = 1.17 for 
ideal gas (y = 5/3) [14], Rs, the temperature T , and the shock velocity v; are written 
as 


Ex, NUS 
r, - so( =) t3) (pc) , (7) 
no 
10 f £51 | 5-3 
T, = 1.5x10 (Sn (K) (8) 
no 
2/5 
= 1.2x 108 x (=) t75 (K), (9) 
no 
R E 1/5 
vV; = dR; =2.1x 10° x (=) RE (cm s7!) : (10) 
dt no 


where E51, no, and t3 are the explosion energy in the unit of 10°! erg, the number 
density of the ambient gas in the unit of cm~’, and the SNR age in the unit of 
10° years. During the adiabatic expansion phase, ~ 70% of the initial explosion 
energy is converted into thermal energy of swept-up matter [15]. 

Remnants in the adiabatic or Sedov expansion phase are typically bright at X-ray 
wavelengths, owing to both the hot ISM gas shocked by the forward shock and the 
hot ejecta shocked by the reverse shock. In Figure 3, we show a canonical example 
of a remnant in this phase: DEM L71, a remnant in the LMC believed to be a few 
thousand years post-explosion [16, 17]. 

High-mass stars typically eject mass as stellar winds before their explosions, 
forming a circumstellar medium (CSM). The winds make a density gradient around 
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Fig. 3 SNR DEM L71, located in the LMC. This X-ray image shows a clear separation between 
the forward shocked material in the ring on the outside and the reverse-shocked material in blue on 
the inside. Credit: NASA/CXC/Rutgers/J.Hughes et al. 


the progenitor star, making the deceleration of the shock slower. A modified simi- 
larity solution for this case is shown in [18]. 

These models introduced here ignored radiative cooling under the assumption 
of an ideal gas. In the case that the system is not truly adiabatic, energy is robbed 
from the shock and v; becomes smaller. Such a case can be caused by efficient 
acceleration of cosmic rays at the shock front [19, for example]. 


3.3 Snowplow phase 


When the shock velocity has decelerated to ~200 km s^, the temperature of the 
gas drops below ~ 10° K and radiative cooling begins to affect the dynamics of 
the shock evolution. Cooling of the gas is a runaway process, as the more it cools, 
the more the cooling rate increases. The density of the shock front becomes highest 
where the timescale of radiative cooling becomes shortest; as a result, a cool, dense 
shell is formed just behind the shock front, whereas the inner gas still has high 
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temperature and high pressure due to the low density. It is during this phase that the 
remnant becomes bright in optical emission from recombining atoms, most notably 
Ha, [S II], [N II], and [O III]. As an example of a remnant in this phase, we show 
an HST image of a portion of the Cygnus Loop in Figure 4. The Cygnus Loop is the 
remains of a SN believed to be 10,000-20,000 yr post explosion. 


Fig. 4 A portion of the Cygnus Loop SNR, seen with HST. Narrowband images of [O III], [S II], 
and Ha, plus V and / continuum bands combine to make this image. Credit: NASA/ESA/Hubble 
Heritage Team. 


The cool shell continues the expansion further collecting ISM gas in a manner 
reminiscent of a snowplow (hence the name "snowplow phase"). At first, the gas 
expands adiabatically, then pV5/ ? becomes a constant. The shock expands with the 
time dependency of R, œ 12/7, As the temperature cools down further, the pressure 
can be ignored and the cool shell expands at a constant radial momentum with a 
time dependency of roughly R; œ yu 


3.4 Dissipation phase 


As the SNR gets even older, the expansion velocity becomes smaller. When the 
shock slows down to a velocity comparable to the sound velocity, the SNR blends 
into the interstellar medium. This is called the “Dissipation phase,” and marks the 
end of the SNR phase. The lifetime of typical SNRs is 106 yrs, during which time 
they will expand to radii of dozens to hundreds of parsecs, filling a volume large 
enough to encompass hundreds to thousands of star systems. 
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4 Types of supernova remnants 


4.1 Shell-type SNRs 


The majority of remnants are classified as "shell-type" SNRs, where the forward 
shock is delineated by emission at wavelengths across the electro-magnetic spec- 
trum. The remnants shown in Figures 1—4 are all shell-type SNRs. The physical 
processes responsible for emitting at various wavelengths differ. In radio regime, 
emission is dominated by synchrotron emission from energetic electrons (acceler- 
ated by the shock wave to GeV energies) spiraling around the turbulent magnetic 
field at the shock front. Historically, most SNRs have been discovered and catego- 
rized based on their radio emission [1]. 

At infrared wavelengths, young SNRs (those in the free expansion or adiabatic 
expansion phases) are generally dominated by thermal emission from warm dust 
grains; see [26] for a review. Dust grains themselves do not generally feel the pas- 
sage of the shock wave, but are suddenly immersed in the hot plasma behind the 
shock, where ion and electron temperatures can be tens of millions of degrees. 
Grains are collisionally heated by these particles, and reach temperatures of 50-200 
K. Dust at these temperatures radiates in the mid-IR, with thermal spectra peaking 
in the 20-100 um range. 

Generally speaking, IR emission from young shell-type SNRs only shows up in 
emission from the shell of the remnant. This emission results from warmed interstel- 
lar grains; searches for dust associated with the ejecta in SNRs has turned up very 
little [27, 28]. As an example, in Figure 5, we show the IR emission from DEM L71, 
the remnant whose X-ray emission is shown in Figure 3. 


Fig. 5 Left: The Spitzer 24 um emission from DEM L71. Middle: The remnant in Hæ. Right: The 
remnant in X-rays. Figure modified from Figure 1 in [27]. 


At optical wavelengths, emission from the shell of a shell-type remnant can be 
seen in one of two ways. If there is at least partially neutral material present ahead 
of the shock (often taken as an indication that the object is the remnant of a Type Ia 
SN; see [29]), then charge exchange can take place in the post-shock environment 
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giving rise to a strong Hæ component [30], as shown in the middle panel of Fig- 
ure 5. The temperatures in these remnants is too hot to allow recombination lines in 
other elements. In older remnants where the temperatures have cooled significantly, 
recombination lines from various elements begin to dominate (see Section 3). 

At X-ray wavelengths, shell-type SNRs may show emission from either or both 
of the forward-shocked ISM/CSM and the reverse-shocked ejecta, depending on the 
evolutionary state of the remnant. The physical processes causing the X-ray emis- 
sion can be either thermal, non-thermal, or both. The thermal emission is generally 
a combination of both Bremmsstrahlung continuum produced when hot ions and 
electrons interact, as well as line emission from highly-ionized atoms of various 
metals, such as O, Fe, Ne, Si, and S. Nonthermal emission in SNRs arises from 
synchrotron emission identical to that seen in radio waves, but from much more en- 
ergetic (~ 10 — 100 TeV) electrons. Because electrons can only be accelerated to 
energies this high very close to the blast wave, synchrotron emission is most often 
seen from the outer edges of the shell in shell-type SNRs [31, 32]. 

Shell-type remnants with sufficient spatial extent can even be identified in 
gamma-rays. For example, [33] used Fermi observations to show that the Cygnus 
Loop is spatially-resolved into a shell in the 0.5-10 GeV energy range. The various 
physical processes involved in gamma-ray emission from SNRs is beyond the scope 
of this chapter; for a thorough review, see [34]. 


Fig. 6 X-ray image of the Crab nebula taken with Chandra. (credit: NASA/CXC/SAO) 
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4.2 Plerion-type SNRs 


Young neutron stars born in core-collapse SNe often form pulsar wind nebulae. 
Pulsar wind nebulae emit bright nonthermal emission from radio to very high energy 
gamma-rays, via synchrotron or inverse Compton processes. These are referred to 
as "plerion-type," or “plerionic” SNRs [35]. 

The typical plerion-type is the Crab nebula. It is the remnant of SN 1054, and 
is now one of the brightest celestial sources at most wavelengths. Figure 6 shows 
the X-ray image of the Crab nebula taken by Chandra. Jets extend out to northwest 
and southeast from the central source (the pulsar), while two torii are seen ema- 
nating outward from the center. The inner radius of the torus is 0.01—0.1 pc, which 
corresponds the shock of the relativistic winds from the pulsar terminated by am- 
bient material [20, 21, 22]. For more details, please see the chapter of this work 
concerning PWNe. 

Some plerion-type SNRs are surrounded by shells of emission resulting from the 
forward shock interacting with the surrounding medium. Such remnants are often 
called “composite type.” One of the typical composite type SNRs is G21.5—0.9 
[23]. Figure 7 is the X-ray image of G21.5—0.9. One can see both the central PWN 
and the outer shell. 

The Crab nebula is a case of a remnant without its shell; it is still an open issue 
why there is no detected shell surrounding the Crab nebula. Recent observations 
by the Hitomi satellite showed that the thermal X-ray emission around the Crab 
nebula is significantly fainter than other young SNRs [24]. The total ejecta mass 
is estimated to be less than ~ 1Mo, implying that the progenitor explosion was 
relatively low energy. One possibility is that the Crab nebula is a remnant of an 
electron capture supernova [36]. More detailed observations for the Crab and more 
samples are needed to understand the origin of plerion-type SNRs. 


4.3 Mixed-Morphology SNRs 


Some relatively evolved SNRs show different morphology in radio and X-ray bands; 
shell-like structure in radio, whereas the centrally-filed in X-rays. They are called 
“Mixed-Morphology SNRs” or MM SNRs [25]. The X-rays are thermal, thus it is 
not due to its central pulsar and/or pulsar wine nebula like plerion cases. 

The mechanism by which mixed-morphology remnants are formed is still un- 
clear. Several hypotheses have been discussed to explain how a remnant with a radio 
shell can appear centrally filled in X-rays. The first possibility is “fossil emission" 
of the remnant [37]; the outer shell cools down due to the expansion and radiative 
cooling, and as a result, the shock temperature becomes too low to emit thermal 
X-rays. Synchrotron emission from the accelerated electrons in the shock disappear 
from the high energy side due to the synchrotron cooling. Thus the shells emitting 
synchrotron disappear within a few thousand years in the X-ray band, whereas they 
remain for more than ten thousand years in the radio band. 
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Fig. 7 X-ray image of the plerion-type SNR, G21.5—0.9 taken with Chandra. (credit: Credit: 
Heather Matheson Samar Safi-Harb (Univ. Manitoba), CXC, NASA) 


Another possibility is due to “evaporated cloud"; the enhanced interior X-ray 
emission could arise from gas evaporated from clouds [38] in the cases that shocks 
of SNRs propagates through a cloudy interstellar medium. Dense and small size 
clouds can pass the shock and later evaporate and are heated, and as a result makes 
thermal X-ray emission inside the shells. 


5 Conclusions 


Supernova remnants are an incredibly diverse class of objects. From simply a purely 
morphological standpoint, their are influenced by a number of factors, including 
mass loss from the progenitor system, the energy and isotropy of the explosion itself, 
and the structure of the surrounding ISM/CSM. The radiation they emit also varies 
from source to source: some remnants are bright in radio waves, some in optical, 
some in X-rays. In general, the remnant's state of dynamical evolution (which can 
be roughly thought of as the product of age and the density being encountered by 
the forward shock) is the most important factor for determining the spectral energy 
distribution for a particular SNR. Other chapters within this work will go into more 
detail on many aspects of SNR study. 


12 Aya Bamba and Brian J. Williams 


References 


[1] Green, D. A. 2019, Journal of Astrophysics and Astronomy, 40, 36. 
doi:10.1007/s12036-019-9601 -6 
[2] Case, G. L. & Bhattacharya, D. 1998, THe Astrophysical J., 504, 761. 
doi:10.1086/306089 
[3] Maggi, P., Haberl, E, Kavanagh, P. J., et al. 2016, A&A, 585, A162. 
doi:10.1051/0004-6361/201526932 
[4] Sasaki, M., Pietsch, W., Haberl, E, et al. 2012, A&A, 544, A144. 
doi:10.1051/0004-6361/201219025 
[5] Garofali, K., Williams, B. F., Plucinsky, P. P., et al. 2017, MNRAS, 472, 308. 
doi:10.1093/mnras/stx 1905 
[6] Reynolds, S. P., et al. 2008, ApJ, 680, 41 
[7] Koyama, K., Petre, R., Gotthelf, E. V., et al. 1995, Nature, 378, 255 
[8] Aharonian, F. A., Akhperjanian, A. G., Aye, K.-M., et al. 2004, Nature, 432, 
75 
[9] Reynolds, S. P. 2008, ARAA, 46, 89 
[10] Hamilton, A. J. S. & Sarazin, C. L. 1984, The Astrophyical J., 287, 282. 
doi:10.1086/162687 
[11] Hamilton, A. J. S. & Sarazin, C. L. 1984, The Astrophysical J., 281, 682. 
doi:10.1086/162145 
[12] Hughes, J. P. 1999, The Astrophysical J., 527, 298. doi:10.1086/308082 
[13] Sedov, L. I. 1959, Similarity and Dimensional Methods in Mechanics, New 
York: Academic Press, 1959 
[14] Landau, L. D. & Lifshitz, E. M. 1959, Course of theoretical physics, Oxford: 
Pergamon Press, 1959 
[15] Chevalier, R. A. 1974, The Astrophysical J., 188, 501. doi:10.1086/152740 
[16] Hughes, J. P., Ghavamian, P., Rakowski, C. E., & Slane, P. O. 2003, ApJ, 582, 
95 
[17] Rakowski, C. E., Ghavamian, P., & Hughes, J. P. 2003, ApJ, 590, 846 
[18] Truelove, J. K. & McKee, C. F. 1999, The Astrophysical J. Suppliment, 120, 
299. doi:10.1086/313176 
[19] Berezhko, E. G., Petukhov, S. I., & Taneev, S. N. 2002, Astronomy Letters, 
28, 632. doi:10.1134/1.1505508 
[20] Rees, M. J. & Gunn, J. E. 1974, MNRAS, 167, 1 
[21] Kennel, C. F. & Coroniti, F. V. 1984, ApJ, 283, 694 
[22] Bamba, A., Mori, K., & Shibata, S. 2010, ApJ, 709, 507 
[23] Slane, P., Chen, Y., Schulz, N. S., et al. 2000, ApJL, 533, L29 
[24] Hitomi Collaboration, Aharonian, F., Akamatsu, H., et al. 2018, PASJ, 70 
[25] Rho, J. & Petre, R. 1998, ApJL, 503, L167. doi:10.1086/311538 
[26] Dwek, E., & Arendt, R. 1992, ARA&A, 30, 11 
[27] Borkowski, K., et al. 2006, ApJ, 642, 141 
[28] Williams, B.J., et al. 2006, ApJ, 652, 33 
[29] Williams, B.J., et al. 2011, ApJ, 741, 96 
[30] Chevalier, R.A., Kirshner, R.P., & Raymond, J.C. 1980, ApJ, 235, 186 


Supernova remnants: Types and evolution 


1] Bamba, A., Yamazaki, R., Ueno, M., et al. 2003, ApJ, 589, 827 

2] Bamba, A., Yamazaki, R., Yoshida, T., et al. 2005, ApJ, 621, 793 

3] Katagiri, H., et al. 2011, ApJ, 741, 44 

4] Reynolds, S.P. 2008, ARA&A, 46, 89 

5] Weiler, K.W. & Panagia, N. 1978, A&A, 70, 419 

6] Miyaji, S., Nomoto, K., Yokoi, K., & Sugimoto, D. 1980, PASJ, 32, 303 
7] Seward, F. D. 1985, Comments on Astrophysics, 11, 15 

8] 


[3 
[3 
[3 
[3 
[3 
[3 
[3 
[38] White, R. L. & Long, K. S. 1991, ApJ, 373, 543 


